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This work reviews and examines two particular issues related with the new technique of electrical 
detection of ferromagnetic resonance (FMR). This powerful technique has been broadly applied for 
studying magnetization and spin dynamics over the past few years. The first issue is the relation 
and distinction between different mechanisms that give rise to a photovoltage via FMR in composite 
magnetic structures, and the second is the proper analysis of the FMR line shape, which remains the 
"Achilles heel" in interpreting experimental results, especially for either studying the spin pumping 
effect or quantifying the spin Hall angles via the electrically detected FMR. 

PACS numbers: 85.75.-d, 75.40.Gb, 76.50.+g, 42.65.-k 



C/3 



o 



> 
cn 

rn 

in 
o 



X 



I. INTRODUCTION 



Electrical detection of ferromagnetic resonance (FMR) 
in ferromagnets (FM) is a powerful new experimen- 
tal tool which has transformed the research on spin 
and magnetization dynamics.---- Over the past few 
years, this technique has generated a great deal of in- 
terest in the communities of magnetism, spintronics, 
and microwave technologies. It has been broadly ap- 
plied for studying diverse material structures, rang- 
ing from ferromagnetic thin films such as Py (permal- 
loy, Ni8oFe2o)^'^'i^'i^, CrOaji^ Fe304fi^ single crystal 
Fe^ GaMnAsfi^ and Lai_a,Sra;Mn03fiS bilayer devices 
such as Py/Pt^i^^i2i^i26 Py/Au^>2i Py/GaAs,22 and 
Y3Fe50i2/Pt^i^ to a variety of magnetic tunneling 
junctions (MTJ) based on magnetic multilayers 4i2r— 
From a technical standpoint, its high sensitivity has made 
it possible to quantitatively determine spin boundary 
conditions^! and to directly measure non-linear magne- 
tization damping2^"3Q,^ the quasiparticle mass for the do- 
main walP^, the phase diagram of the the spin-transfer 
driven dynamics^ and various kinds of parametric spin 
wave excitation^i^iS. Its capability to probe the in- 
terplay of spins, charges, and photons has been utilized 
for studying spin rectificatio n^^'^^ , spin pumping^, spin 
torquei^, and spin Hall effect o^^'^^i^^ , which have led to 
the proposing and realization of novel dynamic spintronic 
devices such as the spin battery^ii^liS^ spin diode^ '^°d^ 
spin dynamo, ^^'^^ and spin demodulator'^''. Very re- 
cently, its ability to detect coherent processes^"— has 
enabled electrical probing of the spin-resonance phase 
and the relative phase of electromagnetic wavesS^, which 
pave new ways for microwave sensing^, non-destructive 
imaging)^ and dielectric spectroscopy22. Such a coher- 
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ent capability is especially exciting as it resembles the 
latest achievement in semiconductor spintronics, where a 
new platform for coherent optical control of spin/charge 
currents has been developed by using nonresonant quan- 
tum interferences 

From the physical standpoint, many different effects 
may generate a time- independent dc voltage in magnetic 
materials via the FMR. Reported mechanisms involve 
spin rectificatioE^'i^, spin pumping'', spin torque^-, 
spin diode^ii^iii, spin Hall2^ and inverse spin Hall 
effects^i^'2ii^. Two major issues stand out here: (1) 
A unified picture clarifying the relations and distinctions 
between such diverse mechanisms has not been estab- 
lished, which leads to increasing controversy and con- 
fusion in interpreting and understanding experimental 
results. A stunning example of this issue is found in 
the very recent studies of the spin Hall effect via elec- 
trically detected FMR, where two similar experiments 
performed on similar devices were interpreted completely 
differently.^'''^^ (2) When more than one mechanism si- 
multaneously plays a role in the FMR generated dc volt- 
age, proper interpretation requires a quantitative anal- 
ysis of the FMR line shape. In our opinion, this has 
remained the "Achilles heel" in recent studies of spin 
pumping and the spin Hall effect which utilize electrically 
detected FMR. The purpose of this article is to address 
these two critical issues with a brief review of the key 
physics of this subject, followed by systematically mea- 
sured experimental data with detailed theoretical analy- 
sis. 

This paper is split into three main sections. First we 
provide a brief review of different mechanisms which may 
generate the photovoltage via the FMR. Then we use 
the dynamic susceptibility obtained from a solution of 
the Landau-Lifshitz-Gilbert equation to derive analyti- 
cal formulae for analyzing the line shape and the sym- 
metry properties of the photovoltage generated through 
spin rectification. Finally we present experimental re- 
sults measured from different samples, at different fre- 
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FIG. 1: (color online). Dynamic response of magnetic struc- 
tures under microwave irradiation: (a) Single thin film layer 
where the spin rectification is due to the magnetic field torque 
as shown in (e). (b) Magnetic bilayer device which has two 
rf currents j and with different spin polarizations. There- 
fore spin rectification is due to both magnetic field and spin 
torques, (c) Magnetic tunneling junction with both j and j^. 
(d) Coordinate system for single ferromagnetic microstrips 
measured in this work under an in-plane applied static mag- 
netic field H. The z'-axis is fixed along the strip and the 
direction of current flow, while the 2-axis is rotated to follow 
the direction of H. (e) Components of magnetic field torque, 
(f) Spin torque in magnetic tunneling junction. 



quencies, and in different experimental configurations, 
showing that the FMR line shape is determined by the 
relative phase of microwaves which is sample and fre- 
quency dependent. 



II. A BRIEF REVIEW OF ELECTRICAL 
DETECTION OF FMR 

Under microwave excitation at angular frequency oj, 
the rf electric (e) and magnetic (h) fields inside a fer- 
romagnetic material can be described as e = eoe~*'^* 
and h = hoe~*("*~*', respectively. Note that in general, 
due to the inevitable losses of microwaves propagating 
inside the ferromagnetic material, there is a phase dif- 
ference $ between the dynamic e and h fields. Such a 
relative phase is determined by the frequency-dependent 
wave impedance of the materials^. As shown in Fig. 
[U the rf e field drives a rf current j = ere, while the rf 
h field exerts a field torque on the magnetization and 



drives it to precess around its equilibrium direction [Fig. 
[IJe)]. Such a magnetization precession is described by 
the non-equilibrium magnetization m = xh. Here tr and 
X are the high-frequency conductivity and Polder ten- 
sor, respectively. Note that due to the resonance na- 
ture of the precession, m lags h by a spin resonance 
phase 9. However, despite the phase of $ and 0, the 
dynamic j and m keep the coherence of their respec- 
tive driving fields, so that the product of any combi- 
nation of their components may generate a time inde- 
pendent signal proportional to {Re{j) ■ Re{m)), where () 
denotes the time average. The amplitude of such a sig- 
nal depends on the phase difference of j and m, which 
can be easily understood from the trigonometric rela- 
tion: cos{ujt) ■ cos(wi - $) = [cos($) -I- cos(2i:ji - $)]/2. 
This is the spin rectification^^ as we highlight in Table 
I. For transport measurements on magnetic structures 
under microwave irradiation, various magnetoresistance 
effects such as anisotropic magnetoresistance (AMR) , gi- 
ant magnetoresistance (GMR) and tunneling magnetore- 
sistance (TMR) make corrections to Ohm's law via their 
corresponding magnetoresistance term a^^i^^ . Such non- 
linear terms typically lead to the product of j and m. 
Spin rectifications induced by such magnetoresistance ef- 
fects are listed in Table I by the terms labeled Vmr- The 
general feature of Vjsm is that its amplitude depends on 
both the relative phase $ and the spin resonance phase 
Q, which leads to a characteristic phase signature of the 
FMR line shape^li^s. 

Similar to the effect of the rf h field torque, a spin 
torque induced by a spin polarized current may also drive 
magnetization precession. For example, in a bilayer [Fig. 
[IJb)] made of a ferromagnetic layer and a nonmagnetic 
layer with spin-orbit coupling2^, in addition to the rf cur- 
rent j flowing in the ferromagnetic layer, the rf e fleld also 
induces a rf charge current flowing in the nonmagnetic 
layer. Via the spin Hall effect in such a nonmagnetic 
layer with spin-orbit coupling, the rf charge current can 
be converted into a spin current j^,, which may flow into 
the ferromagnetic layer and then drive the magnetiza- 
tion precession via the spin torque. Such a spin torque 
induced non-equilibrium magnetization can be described 
by m = Xjjg, where the spin-torque susceptibility tensor 
Xj introduces a spin resonance phase that is differ- 
ent from & in X. Following a similar consideration for 
the magnetoresistance induced spin rectiflcation, a pho- 
tovoltage depending on the spin Hall effect may be gen- 
erated in the ferromagnetic layer. This is the physical 
origin of the spin Hall induced spin rectiflcation effect,— 
which is listed in Table I by the term labeled Vsh- In 
MTJ [Fig. [Ijc)], the spin polarized current can be 
directly generated in the ferromagnetic layer where the 
magnetization is pinned along a different direction than 
that of the free layer. It tunnels into the free layer and 
drives the magnetization precession via the spin torque 
[Fig. [IJf )] . The induced spin rectification signal has been 
measured in spin diodes^ '"'^"'^^ , which is listed in Table I 
by the term labeled Vsd- 



3 



Over the past few years, systematic studies on spin 
rectifications induced by the field (Vmr) and spin torque 
{VsH, Vsd) have been performed, respectively, at the 
University of Manitobai^^LSOl^^-^^^^ and Cornell 

University^i^iiiii^i^i^. It has been found that due to 
the coherent nature of spin rectification, Vmr , Vsh and 
Vsd all depend on the phase difference between j and 
m. However, only the field torque spin rectification 
(Vmr) can be controlled by the relative phase $ of the 
microwaves.— 

In addition to such coherent spin rectification effects, 
it is known that at the interface between a ferromag- 
netic and a nonmagnetic layer, microwave excitation may 
generate a spin polarized current flowing across the in- 
terface via the spin pumping effect^*. This effect has 
been observed in a few striking experiments by measur- 
ing either transmission electron spin resonance^'' or en- 
hanced magnetization damping*^. It involves FMR, ex- 
change coupling and non-equilibrium spin diffusion. In 
our opinion the physical picture of spin pumping was 



best explained in the classical paper of Silsbee et.al. [Ref. 
|46| . which highlighted the key mechanism of dynamic 
exchange coupling between the precessing magnetization 
and the spin polarized current. Such a dynamic cou- 
pling significantly " amplifies" the effect of the rf h field 
in generating non-equilibrium spins. It was later pro- 
posed that the spin current generated via spin pumping 
may also induce a photovoltage, either across the inter- 
face in a spin batterjiMi^^^ or within the nonmagnetic 
layer via the inverse spin Hall effect^ i^"'^^'^^ . Recent ex- 
periments performed on magnetic bilayers^^ have found 
that spin-pumping induced dc voltage (the term Vsp in 
Table I) should be about two orders of magnitude smaller 
than spin Hall induced spin rectification (the term la- 
beled Vsh)- In contrast to phase sensitive coherent spin 
rectification effects, the proposed spin-pumping photo- 
voltage is based on incoherent spin diffusion and FMR 
absorption. Hence, the anticipated FMR line shape is 
symmetric and phase-independent. 



TABLE I: Relation and distinctions between different mech- 
anisms for microwave photovoltages induced by FMR. (For 
simplicity we consider only one matrix element of x and Xj 
which is responsible for the spin rectification, j and rfi de- 
note a corresponding component of the time-dependent cur- 
rent and magnetization, respectively.) 
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Vmr'- Spin Rectification caused by MagnetoResistances^'- 
VsH'- Spin Rectification caused by Spin Hall effect)^ 
Vsd'- Spin Rectification caused by Spm Diode effect •^'iS'ii 
Vsp- Photovoltage caused by Spin Pumping J^^^^"'^^' ' ^^ 
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From the above discussion, it is clear that the line 
shape analysis plays the essential role in distinguishing 
the microwave photovoltage generated by different mech- 
anisms. This issue has been partially addressed by a 
number of theoretical^i^ and experimental work o^i^°d^ 



studying nanostructured MTJs where the photovoltage 
is dominated by the spin torque induced spin rectifica- 
tion. Enlightened by these works and also based on our 
own previous studieai^i^, we discuss in the following the 
critical issue of FMR line shape analysis in microstruc- 
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tured devices, where the field and spin torque induced 
spin rectification may have comparable strength. Our 
theoretical consideration and experimental data demon- 
strate the pivotal role of the relative phase which 
was often under-estimated in previous studies. Via sys- 
tematic studies with different device structures, mea- 
surement configurations and frequency ranges, we find 
that $ has to be calibrated at different microwave fre- 
quencies for each device independently. Hence, our re- 
sults are in strong contradiction with the recent exper- 
iment performed on microstructured magnetic bilayers 
for quantifying the spin Hall angles, where $ was de- 
clared to be zero for all devices at different microwave 
frequencies2Si2i. 



III. FMR LINE SHAPE 



A. The Characteristic Signature 



From Table I, the role of the phase in the FMR line 
shape symmetry can be understood by considering the 
spin rectified voltage V cx {Re{j) ■ Re{rh)). For spin 
rectification induced by the field torque, depending on 
the experimental configuration, at least one matrix com- 
ponent X of the Polder tensor x will drive the FMR; 
whether an on or off-diagonal component is responsi- 
ble for the magnetization precession depends on the 
measurement configuration. Since m = xh, Re{rh) cx 
Re{x) cos{ujt — <i>) + Im{x) sin{ujt — <&). Therefore after 
time averaging a time independent dc voltage is found 
V{^) cx [i?e(x) cos(<i>) - /to(x) sin(<I>)]. It is weh known 
that for diagonal matrix elements, i?e(x) has a disper- 
sive line shape while /m(x) has a symmetric line shape. 
However since the on and off-diagonal susceptibilities dif- 
fer by a phase of 7r/2, if the FMR is driven by an off- 
diagonal susceptibility, the roles are reversed and Re{x) 
has a symmetric line shape while Im{x) has a dispersive 
line shape. 

Based on the simple argument leading to the above 
V{^) expression, one can see that the line shape symme- 
try has a characteristic dependence on the relative phase 
$ between electric and magnetic fields. Thus when mea- 
suring FMR based on the field torque induced spin rec- 
tification effect, it is important to consider the relative 
phase, whereas for a spin pumping measurement which 
measures |nip, or for a spin torque induced spin recti- 
fication which involves the relative phase does not 
influence the experiment. In the next two sections, a 
detailed analysis is given by solving the Landau-Lifshitz- 
Gilbert equation, which leads to analytical formulae de- 
scribing the symmetric and dispersive line shapes for dif- 
ferent measurement configurations. 



B. The Dynamic Susceptibility 

The Landau-Lifshitz-Gilbert equation provides a phe- 
nomenological description of ferromagnetic dynamics 
based on a torque provided by the internal magnetic field 
Hi which acts on the magnetization M, causing it to 



precesi 



,50 



— = -7(MxH.) + - 



MX — ). (1) 



Here 7 is the effective electron gyromagnetic ratio and 
a is the Gilbert damping parameter which can be used 
to determine the FMR line width AiJ, according to 
A_ff ~ auj/^. For the case of microwave induced fer- 
romagnetic resonance Eq. ([ij can be solved by splitting 
the internal field into dc and rf components and tak- 
ing the applied dc field H, along the z-axis. We can 
relate the internal field Hj = Hoi -I- h^e"*"*, to the 
applied field through the demagnetization factors iV^, 
^^Ojz = H ~ N^Mq, hik = hke^'^'' - Nknik, where $fc is 
the relative phase shift between the electric and magnetic 
fields in the k*^ direction and Mq is the dc magnetization 
also along the z-axis. With the magnetization separated 
into dc and rf contributions M — Mq -I- me~"^*, the so- 
lution of Eq. ([IJ yields the dynamic susceptibility tensor 
X which relates the magnetization m to the externally 
applied rf field h 



m = xh 



~^Xxy Xyy ' h 



\Xxx\ ^ 

Xxy\e-'^ 





(2) 



where Q = arctan[Ai//(iJ — Hr)] is the spin resonance 
phasc'^" which describes the phase shift between the re- 
sponse and the driving force in terms of the line width 
AH and the resonance field Hr which are constant for a 
fixed frequency. Q will change from 180° (driving force 
out of phase) to 0° (driving force in phase) around the 
resonance position, in a range on the order of A.H , pass- 
ing through 90° at resonance. This represents the uni- 
versal feature of a resonance; the phase of the dynamic 
response always lags behind the driving force.^^ 

To emphasize the resonant feature of the susceptibil- 
ity tensor elements we define the symmetric Lorentz line 
shape L, and the dispersive line shape D as 



L 



D = 



[H - Hr f + AiJ2 ' 

AiJ(iJ - Hr) 

[H - Hrf + AH^ ■ 



(3) 



Clearly the spin resonance phase can also be written 
in terms of L and as 9 = arctan[A_ff/(iJ — Hr)] = 
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arctan(L/_D) so that L oc sin(0) and D oc cos(O). There- 
fore L and D carry the resonant information of the sus- 
ceptibiUty tensor. 

Using L and D aUows the elements of x to be written 

&S {Xxx T Xxy , Xyy) ~ {D ~\- iL) i^Axx, Axy, Ayy). AxxjA 



xy 



and Ayy are real amplitudes which are related to the 
sample properties 



Ax 



A — _ 
-^xy 



A — 

^yy ~ 



jMojMoNy + {H~ N,Mo)) 
auj{2{H - N,Mo) + Mo{N, + Ny)) ' 
Mo 



a(2(i7 - N,Mo) + Mo{N, + Ny)) ' 
-fMoiMoNx + {H- N,Mo)) 



aLo{2{H - N,Mo) + Mo{N, + Ny)) ' 



(4) 



Since these amplitudes are real all components of x 
include both a dispersive and a Lorentz line shape deter- 
mined solely from the D -\-iL term. However, in a trans- 
mission experiment performed using a resonance cavity 
oc + = L IS measured. This product removes 
the phase dependence carried by L and D and leaves 
only the Lorentz line shape. For the same reason, the 
microwave photovoltage induced by spin pumping (the 
Vsp term in Table I) has a symmetric line shape. 

The susceptibility for the two cases of in-plane and 
perpendicularly applied dc magnetic fields can easily be 
found from Eq. ^ by using the appropriate demagne- 
tization factors. When the lateral dimensions are much 
larger than the thickness, N^; = N2 = and Nj, = 1 
for an in-plane field and N^^ = Nj^ = and = 1 for 
a field applied at a small angle from the perpendicular. 
In this paper, we focus on the in-plane case. The line 
shape analysis for the perpendicular case can be found 
in Ref. |33- In both cases the form of the susceptibility, 
X oc D + iL, describes the ferromagnetic resonance line 
shape where each element of x is the sum of an antisym- 
metric and symmetric Lorentz line shape. As we describe 
in the next section, via the Vmr term of the spin recti- 
fication effect, the symmetry properties of the dynamic 
susceptibility influence the symmetry of the electrically 
detected FMR which can be controlled by tuning the rel- 
ative electromagnetic phase 



C. Spin Rectification Induced by the Field Torque 

The field-torque spin rectification effect results in the 
production of a dc voltage from the non-linear coupling 
of rf electric and magnetic fields. For example, it may 
follow from the generalized Ohm's law^^^ 
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(J • M)M + aRnJ x M, 



(5) 



where a is the conductivity, Ap is the resistivity change 
due to AMR and Rh is the extraordinary Hall coefficient. 

As shown in Fig. [21 we use two coordinate systems 
to describe a long narrow strip under the rotating in- 
plane magnetic field H. The sample coordinate system 
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FIG. 2: (color online). Left panel (a) Coordinate system for 
an in-plane dc H field applied along the z-axis at an angle 9h 
with respect to the z'-axis, with a rf /i-field along the a;'-axis. 
(b) The calculated photovoltage (PV) spectrum at 6h = 45° 
and (c) the calculated amplitude of the PV spectrum at FMR 
as a function of Oh according to Eq. (|9}. Right Panel (d)- 
(f) are the same as (a)-(c), respectively, but with a rf /i-field 
along the y-axis, and calculations are according to Eq. (I10|l . 
In both cases, $ is assumed to be zero. 



(x',y,z') is fixed with the sample length along the z' 
direction and the sample width in the x' direction. The 
measurement coordinate system (x,y^ z) rotates with the 
H direction which is along the z axis. We define Oh as 
the angle between the direction of the strip and the in- 
plane applied static magnetic field (i.e., between the z' 
and z directions). In both coordinate systems, the y axis 
is along the normal of the sample plane. In the case of a 
sample length much larger than the width, the rf current, 
j = Jz'e~'"' flows along the strip direction z' . In this 
geometry the field due to the Hall effect will only be in 
the transverse direction and will not generate a voltage 
along the strip. Taking the time average of the electric 
field integrated along the z' direction, the photovoltage 
is found asi^ii^ 



V : 



AR 

Mb 



{Re{j)-Reimx))sm{29H), 



(6) 



where AR is the resistance change due to the AMR effect 
and the sin(26'i/) term is a result of the AMR effect which 
couples J and M. 

The susceptibility tensor given by Eqs. ^ and ^ 
can be used to write ifix in terms of the rf h field. Since 
Mo and H are both along the z-axis, only the compo- 
nents of h perpendicular to z will contribute to m. How- 
ever, since the rf current flows in the z' direction, to 
calculate the rectified voltage, rhx must be transformed 
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into the {x',y,z') coordinate system by using the rota- 
tion (x,y,z) = (cos(6'//)x' — sin(6'if )z', y, sin(6'jy)x' + 
cos(0/f)z'), which introduces an additional 6h depen- 
dence into the photovoltage. We find that the photo- 
vohage can be written in terms of the symmetric and 
antisymmetric Lorentz line shapes, L and D, as 

A R 

V^^^j,,{AlL + AdD), (7) 

where 

Al = sm{2e H)[-A^^h^> cos{0h) sin{'^^,) 

cos{^y) + A^^h^, sin(6'^f) sin($2/)], 

Ad = sm{2d H)[Axxhx' cos(6'if ) cos{^x') 

sin($y) - Axxhz' sm{9H) cos(<I'z')]) 

(8) 

and 4>^/, $j, and are the relative phases between elec- 
tric and magnetic fields in the x\y and z' directions, 
respectively. 

The amplitudes of the Lorentz and dispersive line 
shape contributions show a complex dependence on the 
relative phases for the x' , y and z' directions and in gen- 
eral both line shapes will be present. However, depending 
on the experimental conditions, this dependence may be 
simplified. For instance when h^' is the dominate driving 
field as shown in Fig. HJa), we may take hy — hz' ~ 
and = ^, which results in 

Ai? 

V = --—;-jz'Axxhx' cos{9h) sm{29H) 

[Lsin($) - i:'cos($)] . (9) 

From Eq. ^ we see that the photovoltage line shape 
changes from purely symmetric to purely antisymmetric 
in 90° intervals of being purely antisymmetric when 
$ = n X 180° and purely symmetric when $ = (2n + 1) x 
90°, n = 0,± 1,± 2.... 

As shown in Fig. Efb) and (c), the photovoltage in 
Eq. ^ also shows symmetries depending on the static 
field direction 9h- Since H -H corresponds to 9h 
9„ + 180°, V{H) = -V{-H). Furthermore at 9h = 
n X 90°, n = 0, ± 1, ± 2 . . . the voltage will be zero. 

Similarly when hy dominates as shown in Fig. (^Jc) , we 
take hx' = hz' ~ Q and ^y = $ which results in a voltage 

Ai? 

^ " " 2A/o-^^"^^^'*^ sin(26lH) 

[Lcos($) + i:isin($)] . (10) 

The symmetry properties are now such that the line 
shape is purely symmetric when ^ — n x 180° and 
purely antisymmetric when $ = (2n + 1) x 90°, n ~ 

0, ± 1,± 2 Also the photovoltage determined by 

Eq. ()10|1 is now symmetric with respect to H under 
0H en + 180° so that V{H) = V{-H) as shown in Fig. 
[5](e). Therefore, experimentally the different symmetry 



of the FMR at H and —H can be used as an indication 
of which component of the h field is dominant. 

Both Eq. © and Eq. pUj) demonstrate that a change 
in the relative electromagnetic phase is expected to result 
in a change in the line shape of the electrically detected 
FMR. It is worth noting that when the relative phase 
<i> = 0, the line shape is purely antisymmetric for FMR 
driven by hx' and purely symmetric for FMR driven by 
hy as illustrated in Fig. [Hb) and[2Ie), respectively. In 
the general case when rhx is driven by multiple h compo- 
nents, Eq. ([7]) must be used in combination with angular 
{9}i) dependent measurements in order to distinguish dif- 
ferent contributions. 



D. The Physics of $ 

It is clear therefore that for field torque induced spin 
rectification, the relative phase $ between the microwave 
electric and magnetic fields plays the pivotal role in the 
FMR line shape. Note that $ is a material and frequency 
dependent property which is related to the losses in the 
system4ii^i^ When a plane electromagnetic wave prop- 
agates through free space the electric and magnetic fields 
are in phase and orthogonal to each other. However 
when the same electromagnetic wave travels through a 
dispersive medium where the wave vector is complex, the 
imaginary contribution can create a phase shift between 
electric and magnetic fields. The most well known ex- 
ample is that of a plane electromagnetic wave moving in 
a conductor*^ where Faraday's law gives a simple rela- 
tion between electric and magnetic fields, w^h = k x e. 
Therefore the complex part of the wave vector k will in- 
duce a phase shift between electric and magnetic fields. 
Although the field will exponentially decay inside a con- 
ductor, it will still penetrate a distance on the order of 
the skin depth, and in a perfect conductor the conduc- 
tivity, which produces an imaginary dielectric constant, 
will result in a phase shift of 45° between the electric and 
magnetic fields 4i 

In a complex system such as an experimental set up in- 
volving waveguides, coaxial cables, bonding wires and a 
sample holder, which are required for electrical FMR de- 
tection, the relative phase cannot be simply calculated. 
Nevertheless losses in the system which can be charac- 
terized in a variety of ways, such as through the wave 
impedance )^^'^^ will lead to a phase shift between elec- 
tric and magnetic fields which will influence the FMR 
line shape. 

Although the physics of $ is in principle contained in 
Maxwell's equations, due to the lack of technical tools 
for simultaneously and coherently probing both e and h 
fields, the effect of the relative phase had often been ig- 
nored until the recent development of spintronic Michel- 
son interferometry^l. In the following we provide sys- 
tematically measured data showing the influence of the 
relative phase $ on the line shape of FMR which is driven 
by different h field components. 



7 



IV. EXPERIMENTAL LINE SHAPE 
MEASUREMENTS 

A. hy Dominant FMR 

In order to use the hy field to drive FMR a first gen- 
eration spin dynamo was used where a Cu/Cr coplanar 
waveguide (CPW) was fabricated beside a Py microstrip 
with dimension 300 /im x 20 fim x 50 nm on a Si02/Si 
substrate as shown in Fig. [HJa). A microwave current 
is directly injected into the CPW and flows in the z' di- 
rection inducing a current in the Py strip also along the 
2;'-axis. In this geometry the dominant rf h field in the 
Py will be the Oersted field in the y direction produced 
according to Ampere's Law. This field will induce FMR 
precession with the same cone angle independent of the 
static H orientation. 

The AMR resistance depends on the orientation of the 
magnetization relative to the current and follows the rela- 
tion R{H) = R{0)~ARsin^{eM), where 0m (not shown) 
is the angle between the magnetization and the current 
direction. For Py the AMR effect, which is responsible for 
the spin rectification, is observed to produce a resistance 
change of AR/R{0) ~ 0.4 %. When H is applied along 
the a;'-axis, i.e., the in-plane hard axis, the magnetization 
M tends to align toward the static field H and the angle 
9m is determined by sin(6'M) = H/Ha for H < Ha, 
where Ha = N^'Mq is the in-plane shape anisotropy 
field. The measured data (symbols) shown in Fig. (JJc) is 
fit (solid curve) according to R{H) = R{0)-ARsm'^{eM) 
with i?(0) = 112.66 n, AR = 0.47 n, hqHa = 4.0 mT, 
and iV^- =0.004. 

Fig. ^d) shows that the line shape at 9h — 120° 
and w/27r ~ 5 GHz is almost purely dispersive, in- 
dicating that at this frequency $ ^ 90° according 
to Eq. ((TU)) . The 9h dependence of Hr is shown in 
Fig. [3fe) and can be well fit by the function lo = 

j^{\Hr\ + Ha cos{29HmHr\ + Mo - Ha{1 + sin" {9h)] 
by taking the shape anisotropy field Ha along the x'-axis 
into account)^ As expected the amplitude of these os- 
cillations is hqHa = 4.0 mT. The frequency dependence 
of Hr at 9h = 45° is shown in Fig. |3l^f) and is fit using 
w = j^\Hr\{\Hr\ +Ma) with 7/27r = 29.0 ^ioGHz/T 
and fioM=1.0 T. 

By systematically measuring the line shape as a func- 
tion of the microwave frequency, we observe the interest- 
ing results of Fig. S) The FMR line shape is observed 
to change from almost purely dispersive at a;/27r = 5 
GHz to almost purely symmetric at a;/27r = 5.56 GHz. 
As discussed before, the line shape may be affected by 
the h orientation, i.e., different h vector components will 
affect the line shape differently. Hence, if changing the 
microwave frequency changes the dominant driving field, 
the line shape may change. To rule out such a possibil- 
ity an angular dependent experiment was performed to 
measure the line shape at different 9jj for each frequency 
uj. The results are plotted on the right panel of Fig. 0] 
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FIG. 3: (color online), (a) Schematic diagram of the first 
generation spin dynamo where the Py strip is located be- 
side the CPW. The dominate magnetic field in the Py is 
the Oersted field in the y direction due to the current in 
the CPW. (b) Micrograph of the device, (c) Magnetoresis- 
tance at 9h = 90°. AMR is seen to be ~ 0.4%. Arrows 
denote the anisotropic field, ^oHa = 4.0 mT. Open circles 
are experimental data and solid curve is the fitting result 
using R(Q) = 112.66 fi, Ai? = 0.47 Q.,Ha = 4.0 mT. (d) 
Electrically detected FMR at Oh = 120° and w/27r = 5 GHz 
showing an almost purely dispersive line shape ($ ~ 90°). 
Fit is according to Eq. (|10p with hqAH = 3.6 mT, noHr 
= 32.2 mT. (e) Oscillating Hr dependence on the static field 
direction 6h with amplitude 2Ha. (f) Dependence of FMR 
frequency on the resonant field Hr at 9h ~ 45° . Open circles 
are experimental data and the solid line is the fit according 
to = -y^lHrlilHrl+Mo). 



which shows the sinusoidal curves for the Lorentz, A]^, 
and dispersive, A]j, amplitudes (dashed and solid curves 
respectively) as a function of the static field angle 9h- 
Both the Lorentz and dispersive amplitudes are found to 
follow a sin(26'^r) dependence on the field angle in agree- 
ment with Eq. (jlOp indicating that the magnetization 
precession is indeed dominantly driven by the hy field. 
Therefore the line shape change indicates that the rela- 
tive phase is frequency dependent. As shown in Fig. 
[5fa), at uj/2tt ~ 5 GHz the amplitude oi A]j is approx- 
imately one order of magnitude larger than A^, while 
at a;/27r = 5.56 GHz Ad is one order of magnitude less 
than Al. Such a large change in Al/Ad shows that in 
a microwave frequency range as narrow as 0.6 GHz, the 
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FIG. 5: (color online), (a) The Ad/Al ratio as a function of 
Lo/2n showing the line shape change from dispersive at 5 GHz 
(left inset) to Lorentz at 5.6 GHz (right inset) with a step size 
of 0.01 GHz. (b) $ dependence on uj/2n over same frequency 
interval showing the same dependence as Ad/Al- 



propagation in a waveguide where the field distribution 
i.e. the waveguide modes, are known to depend strongly 
on boundary conditions and frequency)^ Despite the 
complex wave properties, the key message of our results is 
clear and consistent with the consideration of the physics 
of the relative phase: it shows that in order to properly 
analyze the FMR line shape, $ has to be determined for 
each frequency independently. 



FIG. 4: (color online). Data shown for a first generation spin 
dynamo. FMR spectra at 9h = 120° for several frequencies 
from 5.0 to 5.56 GHz with corresponding Lorentz and dis- 
persive amplitudes as a function of 9h. Circles and squares 
indicate the Lorentz and dispersive amplitudes of Eq. (|10p re- 
spectively and show a sm{26H) dependence as expected. Solid 
and dashed curves are sin(2S£r) functions. 



relative phase <& can change by 90°. Fig. [5{b) shows $ 
determined by using Eq. (jlOp . which smoothly changes 
with microwave frequency except for a feature near 5.18 
GHz, which is possibly caused by a resonant waveguide 
mode at this frequency. 

Such a large change of <i> within a very narrow range 
of microwave frequency indicates the complexity of wave 
physics. Note that microwaves at ^ 5 GHz have wave- 
lengths on the order of a few centimeters which are much 
larger than the sub-millimeter sample dimensions. Con- 
sequently the microwave propagation depends strongly 
on the boundary conditions of Maxwell's equations which 
physically include the bonding wire, chip carrier, as well 
as the sample holder. This is similar to the microwave 



B. hx' Dominant FMR 

In order to drive the FMR using the rf field in the x' 
direction, h^' , a second generation spin dynamo was fab- 
ricated with the Py strip underneath the CPW as shown 
in Fig. ini In this case the 300 fim x 7 fim x 100 nm Py 
strip is underneath the Cu/Cr coplanar waveguide which 
is fabricated on a Si02/Si substrate. Again a microwave 
current is directly injected into the CPW and induces a 
current in the z' direction in the Py strip. The dominant 
rf field in the Py is still the Oersted field, but due to the 
new geometry it is in the x' direction. 

Due to the smaller width and larger thickness, the de- 
magnetization factor, Nx' = 0.008 is twice that in the 
first generation sample. This corresponds to (IqHa = 
8.0 mT as indicated by the broader AMR curve in Fig. 
^c). This value is further confirmed by the Hr vs 6h 
plot shown in Fig. ^e). Fig. ^i) shows the frequency 
dependence of Hr for FMR (circles) and for the first per- 
pendicular standing spin wave resonance (SWR) (trian- 
gles) measured at 9h — 45°. The frequency dependence 

of Hr follows UJ = ^y/{\Hr\ + H^,j,){\Hr\ + A/q -1" ffg^) 

where Hex is the exchange field. In Fig. IHJf ) the standing 
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FIG. 6: (color online), (a) Schematic diagram of the second 
generation spin dynamo where the Py strip is located under- 
neath the CPW. In this case the dominant magnetic field in 
the Py is the Oersted field in the x' direction due to the field in 
the CPW. (b) Micrograph of the Py CPW device, (c) Magne- 
toresistance at 9h = 90°. AMR is seen to be ~ 0.5%. Arrows 
denote the anisotropic field, ^iqHa = 8.0 mT. Open circles are 
experimental data and solid curve is the fitting result using 
R{0) = 121.53 n and AR = 0.66 fi. (d) Electrically detected 
FMR at Oh = 120° and uj/2tt = 8 GHz showing a nearly 
symmetric Lorentz line shape. Fit is according to Eq. (|10|) 
with noAH = 6.0 mT , noHr = 76.5 mT and $ = -102°. 
(e) Oscillating Hr dependence on the static field direction 6h 
with amplitude 2Ha- (f) Dependence of FMR frequency on 
the resonant field Hr at 9h ~ 45°. Solid circles show the 
FMR frequency dependence while the solid triangles are the 
standing S WR frequency d ependence. The solid line is a fit 
to a; = -y^lHrlilHrl+Mo). 



SWR is fit using j/2tt = 29.0 fioGUz/T , noHex = 30 
niT and /UqA/o = 1.0 T. 

Similar to the results presented in the previous section, 
the line shape of FMR measured on the second genera- 
tion sample is also found to be frequency dependent (not 
shown). Hence, $ is found to be non-zero in the general 
case. For example, at uj/2tt — 8 GHz, the line shape is 
found to be nearly symmetric, as shown in Fig. IHJd) for 
the FMR measured at Oh = 120°, which indicates $ is 
close to —90° at this frequency. Note that our result is 
in direct contrast with the recent study of Ref. [l^ and 
[2ll . where experiments were measured in the same con- 




HqH (mT) Cone angle (degree) 

FIG. 7: (color online). Data shown for a second generation 
spin dynamo, (a) FMR line shape at fixed frequency, ijj/2n — 
8 GHz for several Oh from 90° to 180° in steps of 10° . Open 
circles are experimental data and solid lines are fits using Eq. 
(|9]) with <I> — —102° fixed, (b) Ad and Al shown in squares 
and circles respectively as a function of Oh- Fitting curves are 
sin(26'£r) cos{Oh) functions, (c) AH for several values of Oh 
showing an oscillation with Oh- (d) Non-linear dependence of 
line width AH on the cone angle. Dashed line is the expected 
linear Gilbert damping whereas the data follows the quadratic 
dependence shown by the solid line. 

figuration and where it was suggested that $ = 0° for all 
samples at all frequencies. 

While the line shape and hence the relative phase is 
found to be frequency dependent, $ is expected to be 
independent of the static field direction 9h - This is con- 
firmed in Fig. [Tja) which shows the line shape measured 
at several values of 9h in 10° increments. The data can 
be fit well using Eq. (|9]) with a constant $ = —102° for 
all 6h- It confirms that the FMR is driven by a single h 
component, in this case the hx' field, and that $ does not 
depend on 6^- In Fig. [T^b) the 6^ dependence of 
and Ao (sohd/circles and dashed/squares respectively) 
is shown. The circles and squares are experimental data 
while the solid and dashed lines are fitting results using 
a sin(26'if ) cos{9h) function according to Eq. It pro- 
vides further proof that the hx' field is responsible for 
driving the FMR in this sample. 

While the results from both the 1st and 2nd genera- 
tion spin dynamos show consistently that <i> is sample and 
frequency dependent, the 2nd generation spin dynamos 
exhibit special features in comparison with the 1st gen- 
eration spin dynamos: the reduced separation between 
the Py strip and CPW enhances the h.^' field so that the 
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line width /S.H is enhanced by non-hnear magnetization 
daniping25-^.iS^, which depends on the cone angle 9 of the 
precession via the relation 9 ^ h^' cos{9 h) / [9) . As 
shown in Fig. El^c), Ai? is found to oscillate between 4.0 
and 9.0 mT as Oh changes. At 6*// = 0°, 6* /AH and 
the cone angle is at its largest (about 4°). As 9h increases 
from 0° and moves toward 90°, 9 decreases so that the 
non-linear damping contribution to AH decreases. Us- 
ing the cone angle calculated from Fig. [Tjc), we plot 
in Fig. ^d) AH{9) as a function of the cone angle. It 
shows that AH has a quadratic dependence on the pre- 
cession cone angle, which is in agreement with our previ- 
ous study in the perpendicular H-field configuration2Si2^. 
We note that for cone angles above only a few degrees, 
the non-linear damping already dominates the contribu- 
tion to AH . Again, this is in direct contrast with the 
result of Mosendz et al.,—:^, where AH was found to 
be constant by varying 9h, indicating no influence from 
non-linear damping, but the cone angle 9 was estimated 
to be as high as 15° based on the line shape analysis 
assuming relative phase $ = 0. 



C. Arbitrary h Vector 

Next we consider the most general case which is de- 
scribed by Eq. ([T]) where all components of h may con- 
tribute to the FMR line shape. The sample used here is a 
single Py strip where a waveguide with a horn antennae 
provided both the electric and magnetic driving fields. 
The sample chip is mounted near the centre, at the end 
of a rectangular waveguide and the Py strip is directed 
along the short axis of the waveguide. 

In a waveguide, the electromagnetic fields are well 
known and in general three components, h^', hy and hz' 
exist iS Figure ^a) shows both the FMR and perpen- 
dicular standing SWR at 9h = 45°. Indeed both the 
amplitude and the line shape are different for the two 
FMR peaks located at H and —H, which indicates the 
existence of multiple h field components and Eq. ([7]) and 
Eq. (|5]) are needed to separate them. 

This separation is done using the Lorentz and disper- 
sive amplitudes determined from a fit to the FMR which 
are plotted as a function of 9h in Fig. ^h) and (c) for 
uj/2tt = 12 and 11.2 GHz, respectively. A fit using Eq. 

allows a separation of the contributions from each 
of the hx' , hy and hz' fields based on the their different 
contributions to the 9h dependence of the line shape. 

The results of the fit have been tabulated in Table HIl 
where j/27r = 28.0 ^ioGHz/T, ^qMo = 0.97 T and noHr 
— 152 mT were used. The amplitudes of the different h 
field components have been normalized with respect to 
the hx' component. At both 11.2 and 12 GHz the h^' 
field is much larger than hy or hz' , which is expected 
based on the wave propagation in a horn antennae. 

In changing from 11.2 to 12 GHz the relative phase for 
each component is seen to change. Therefore even in the 
case of a complex line shape produced by multiple h field 
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FIG. 8: (color online). Data shown for a single Py strip with 
precession driven by horn antennae field. The strip dimen- 
sions are 3 mm x 50 ^m x 45 nm. (a) Spectra showing dis- 
tinct resonances due to FMR and SWR at uj/2tv — 12 GHz. 
(b) Separated Lorentz and dispersive line shapes (circles and 
squares respectively) as a function of 9h from a fit to Eq. (O 
at uj/2-K = 12 GHz and (c) uj/2n = 11.2 GHz. 



TABLE II: Angular separation of h field components for 12 
and 11.2 GHz. 





12 GHz 


11.2 GHz 


\ha:'\ 


1 


1 


\hy\ 


0.02 


0.14 


\h,'\ 


0.19 


0.37 




-23° 


50° 




40° 


-30° 




-33° 


82° 



components, by separating the individual contributions 
of the rf magnetic field via angular dependence measure- 
ments, the relative phase $ of each field component is 
found to be frequency dependent. 



D. Additional Influences on "1> 

In addition to the frequency and sample dependencies, 
the relative phase $ may also depend on the lead config- 
uration and wiring conditions of a particular device, as 
we have mentioned in Section A. Here, we address such 
additional influences by using the first generation spin 
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FIG. 9: (color online), (a) FMR observed in a first generation 
spin dynamo. Inset shows the device structure with two Py 
strips labeled SI and S2. (b) FMR for Py thickness d = 100 
nm for both SI and S2. In SI $ = -11°, while in S2 the line 
shape is slightly more asymmetric and $ = 22°. (c) For d = 
60 nm the relative phase is $ = -29° for SI and $ = 27° for 
S2. 



dynamosi^ shown in the inset of Fig. IH^a). Two spin 
dynamos with the same lateral dimensions but different 
Py thickness d are studied. Each spin dynamo involves 
two identical Py strips denoted by SI and S2, one in each 
center of the G-S strips of the CPW, which are placed 
symmetrically with respect to the S strip. The current 
and rf h field are induced in the Py via a microwave 
current directly injected into the CPW. Similar to the 
sample discussed in Section A, hy is the dominant field 
which drives the FMR. 

As shown in Fig. EJa), FMR measured at uj/2'k = 5 
GHz on the sample SI with d = 100 nm shows a nearly 
symmetric Lorentz line shape and a field symmetry of 
V{H) = V{-H). From the FMR line shape fitting, $ = 
-11° is found. Interestingly, as shown in Fig. Mjo), the 
FMR of the sample S2 of the same spin dynamo mea- 
sured under the same experimental conditions shows a 
different line shape, from which a different $ — 22° is 
found. We can further compare measured on the other 
spin dynamo with a different Py thickness of c? = 60 nm, 
also at w/27r = 5 GHz. Here for SI, $ = -29° while for 
S2, $ = 27°. Again, the relative phase is found differ- 
ent for SI and S2. These results demonstrate that due 
to additional influences such as a different lead configu- 
ration and wiring conditions, even for samples with the 
same lateral dimensions, $ in each device is not neces- 
sarily the same. It demonstrates clearly that the relative 
phase <i> can not be simply determined by analyzing the 



FMR line shape measured on a reference device. 



E. Closing Remarks 

The experimental data presented above show that re- 
gardless of the FMR driving field configuration, the rel- 
ative phase between the rf electric and magnetic field is 
sample and frequency dependent and non-zero. This non- 
zero phase results in both symmetric and antisymmetric 
Lorentz line shapes in the FMR detected via field-torque 
induced spin rectification. The $ dependence of the line 
shape symmetry changes based on which component of 
the rf h field is responsible for driving the FMR pre- 
cession. For instance a purely antisymmetric line shape 
could correspond to $ = 0° if the FMR is driven by hx', 
or to $ = 90° if the FMR is driven by hy, therefore the 
line shape itself cannot be used to determine $ directly. 
To separate the h field components an angular {9h) de- 
pendence measurement is necessary, which allows both 
h as well as the phase to be determined. Using such a 
measurement $ has been observed to change from 0° to 
90° in a narrow frequency range (0.6 GHz) resulting in 
a change from an antisymmetric to symmetric line shape 
demonstrating the large effect the relative phase has on 
the FMR line shape. Furthermore $ is not identical even 
in samples with the same geometrical size. Therefore in 
our opinion $ cannot be simply determined from a ref- 
erence sample but should be calibrated for each sample, 
at each frequency and for each measurement cycle. 



V. CONCLUSION 

Spin rectifications caused by the coupling between cur- 
rent and magnetization in a ferromagnetic microstrip 
provide a powerful tool for the study of spin dynamics. 
In order to distinguish different mechanisms which en- 
able the electrical detection of FMR via microwave pho- 
tovoltages, it is essential to properly analyze the FMR 
line shape. For spin rectification caused by a microwave 
field torque, due to the coherent nature of this coupling, 
the resulting dc voltage depends strongly on the relative 
phase between the rf electric and magnetic fields used to 
drive the current and magnetization, respectively. There- 
fore not only does electrical FMR detection provide a 
route to study the relative phase, but it also necessi- 
tates calibrating the relative phase prior to performing 
electrically detected FMR experiments. Based on a sys- 
tematic study of the electrically detected FMR, the line 
shape is observed to depend strongly on the microwave 
frequency, driving field configuration, sample structure 
and even wiring conditions. It is in general a combi- 
nation of Lorentz and dispersive contributions. These 
effects have been quantitatively explained by accounting 
for the relative phase shift between electric and magnetic 
fields. Analytical formula have been established to an- 
alyze the FMR line shape. Our results imply that for 



electrically detected FMR which involves both spin Hall 
and spin rectification effects, the pivotal relative phase 
must be calibrated independently in order to properly 
analyze the FMR line shape and quantify the spin Hall 
angle. This cannot be done by using a reference sample 
but could be achieved through such techniques as spin- 
tronic Michelson interferometry^. 
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